Infinite Series 



Infinite series are among the most powerful and useful tools that you've encountered in your introductory 
calculus course. It's easy to get the impression that they are simply a clever exercise in manipulating 
limits and in studying convergence, but they are among the majors tools used in analyzing differential 
equations, in developing methods of numerical analysis, in defining new functions, in estimating the 
behavior of functions, and more. 

2.1 The Basics 

There are a handful of infinite series that you should memorize and should know just as well as you do 
the multiplication table. The first of these is the geometric series, 

oo ^ 

l + x + x"^ + + x'^ -\ = y^x^ = - for Ixl < 1. (2.1) 

It's very easy derive because in this case you can sum the finite form of the series and then take a limit. 
Write the series out to the term x^ and multiply it by {1 — x). 

{1 + x + x^ + x^ + ■■■ + x^){l -x) = 

{l + x + x'^ + x^ + --- + x^)-{x + x^ + x^ + x'^ + --- + x^+^) = 1 - x^+^ (2.2) 

If < 1 then as A*" — )■ oo this last term, x^~^^, goes to zero and you have the answer. If x is outside 
this domain the terms of the infinite series don't even go to zero, so there's no chance for the series to 
converge to anything. 

The finite sum up to is useful on its own. For example it's what you use to compute the 
payments on a loan that's been made at some specified interest rate. You use it to find the pattern of 
light from a diffraction grating. 

yx" = ^^^ (2.3) 

Some other common series that you need to know are power series for elementary functions: 
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\n{l + x) = x-^ + ^-... =E(-1)'''t (N<1) (2-4) 
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Of course, even better than memorizing them is to understand their derivations so well that you 
can derive them as fast as you can write them down. For example, the cosine is the derivative of the 
sine, so if you know the latter series all you have to do is to differentiate it term by term to get the 
cosine series. The logarithm of (1 + x) is an integral of 1/(1 + x) so you can get its series from that 
of the geometric series. The geometric series is a special case of the binomial series for a = —1, but 
it's easier to remember the simple case separately. You can express all of them as special cases of the 
general Taylor series. 

What is the sine of 0.1 radians? Just use the series for the sine and you have the answer, 0.1, or 
to more accuracy, 0.1 - 0.001/6 = 0.099833 

What is the square root of 1.1? y/lJ = (1 + .1)^/^ = 1 + 1 ■ 0.1 = 1.05 

What is 1/1.9? 1/(2 - .1) = 1/[2(1 - .05)] = ^(1 + .05) = .5 + .025 = .525 from the first 
terms of the geometric series. 

What is v^l024? ^1024 = ^1000 + 24 = ^1000(1 + 24/1000) = 
10(1 + 24/1000)1/3 = 10(1 + 8/1000) = 10.08 

As you see from the last two examples you have to cast the problem into a form fitting the 
expansion that you know. When you want to use the binomial series, rearrange and factor your expression 
so that you have 

(l + something small)" 

2.2 Deriving Taylor Series 

How do you derive these series? The simplest way to get any of them is to assume that such a series 
exists and then to deduce its coefficients in sequence. Take the sine for example, assume that you can 
write 

s'mx = A + Bx + Cx"^ + Dx-^ + Ex'^ H 

Evaluate this at x = to get 

sinO = = A + fiO + CO^ + Z^O^ + EQ-^ + --- = A 

so the first term, A = 0. Now differentiate the series, getting 

cos X = 5 + 2Cx + 3Dx^ + AEx'^ H 

Again set X = and all the terms on the right except the first one vanish. 

cos = 1 = + 2C0 + 3D02 + 4^0^ + ... = B 

Keep repeating this process, evaluating in turn all the coefficients of the assumed series. 



sinx 


= A + Bx + Cx^ + Dx^ + Ex^ H 


sinO = 


= 


A 


cosx 


= B + 2Cx + 3Dx^ + 4Ex3 + • • • 


cosO = 


1 = 


B 


sinx 


= 2C + 6Dx + 12Ex^ + ■■■ 


— sin = 


= 


2C 


cosx 


= 6D + 2AEx + 60Fx^ H 


— cos = 


-1 


= 6D 


sinx 


= 24E + 120Fx + • • • 


sinO = 


= 


2AE 


cosx 


= 120F + • • • 


cosO = 


1 = 


120F 



This shows the terms of the series for the sine as in Eq. (2.4). 

Does this show that the series converges? If it converges does it show that it converges to the 
sine? No to both. Each statement requires more work, and I'll leave the second one to advanced 
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calculus books. Even better, when you understand the subject of complex variables, these questions 
about series become much easier to understand. 

The generalization to any function is obvious. You match the coefficients in the assumed expan- 
sion, and get 

fix) = /(o) + xfio) + l^r (0) + |[r (0) + ^r"(o) + • • • 

You don't have to do the expansion about the point zero. Do it about another point instead. 

fit) = fito) + it- to)/'(to) + ^^^r (to) + • • • (2.5) 
What good are infinite series? 

This is sometimes the way that a new function is introduced and developed, typically by determining a 
series solution to a new differential equation. (Chapter 4) 
This is a tool for the numerical evaluation of functions. 

This is an essential tool to understand and invent numerical algorithms for integration, differentiation, 
interpolation, and many other common numerical methods. (Chapter 11) 

To understand the behavior of complex- valued functions of a complex variable you will need to under- 
stand these series for the case that the variable is a complex number. (Chapter 14) 

All the series that I've written above are power series (Taylor series), but there are many other 
possibilities. 

C(^) = E;^ (2-6) 
1 

-'=y + ^l:(-l)"^<=«(^) (-L<^<L) (2^7) 

The first is a Dirichlet series defining the Riemann zeta function, a function that appears in statistical 
mechanics among other places. 

The second is an example of a Fourier series. See chapter five for more of these. 
Still another type of series is the Frobenius series, useful in solving differential equations: its form is 
'^kCtj^x^'^^ . The number s need not be either positive or an integer. Chapter four has many examples 
of this form. 

There are a few technical details about infinite series that you have to go through. In introductory 
calculus courses there can be a tendency to let these few details overwhelm the subject so that you are 
left with the impression that that's all there is, not realizing that this stuff is useful. Still, you do need 
to understand it.* 

2.3 Convergence 

Does an infinite series converge? Does the limit as — > oo of the sum, Uj^, exist? There are a 
few common and useful ways to answer this. The first and really the foundation for the others is the 
comparison test. 

Let U/j and f be sequences of real numbers, positive at least after some value of k. Also assume 
that for all k greater than some finite value, < Vp;. Also assume that the sum, J2k ^fc ^^^^ converge. 
The other sum, Ylik'^k t^^" converges too. This is almost obvious, but it's worth the little effort that 
a proof takes. 



* For animations showing how fast some of these power series converge, check out 
www.physics.miami.edu/nearing/mathmethods/power-animations.html 
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The required observation is that an increasing sequence of real numbers, bounded above, has a 

limit. 

After some point, k = M, all the U). and Vj^ are positive and Uj. < Vj.. The sum an = X^M'^fe 
then forms an increasing sequence of real numbers, so by assumption this has a limit (the series 
converges). The sum bn = Yllvi'^k ^'^ increasing sequence of real numbers also. Because U]^ < Vj, 
you immediately have bn < an for all n. 

bn<an< lim On 

n->oo 

this simply says that the increasing sequence bn has an upper bound, so it has a limit and the theorem 

is proved. 

Ratio Test 

To apply this comparison test you need a stable of known convergent series. One that you do have is 
the geometric series, for \x\ < 1. Let this x^ be the Vf^ of the comparison test. Assume at least 

after some point k = K that all the Uj^ > 0. 
Also that Mfc+i < XU].. 

Then Uj{^2 ^ ^'^K+i u^+i < xuj{ gives ^ x'^uj{ 

You see the immediate extension is 

As long as a; < 1 this is precisely set up for the comparison test using X^^^i^a;" as the series that 
dominates the Yln''^n- This test, the ratio test is more commonly stated for positive Uj^ as 



If for large k, <x<l 



then the series 



converges 



(2.8) 



This is one of the more commonly used convergence tests, not because it's the best, but because it's 
simple and it works a lot of the time. 

Integral Test 

The integral test is another way to check for convergence or divergence. If / is a decreasing posi- 
tive function and you want to determine the convergence of J^nfi''^)' "-^"^ integral 
dx f{x) and check it for convergence. The series and the integral converge or diverge together. 




From the graph you see that the function / lies between the tops of the upper and the lower 
rectangles. The area under the curve of / between n and n + 1 lies between the areas of the two 
rectangles. That's the reason for the assumption that / is decreasing and positive. 



/(n)-l> r^'dxf{x)>f{n + l)-l 

Jn 
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Add these inequalities from n = k to n = oo and you get 

/■fc+l /•k+2 POO 

f{k) + f{k + l) + ...> +/ +•••=/ dxfix) 

Jk Jk+1 Jk 

POO 

>f{k + l) + f{k + 2) + ...> / dxf{x)>f--- (2.9) 

Jk+l 

The only difference between the infinite series on the left and on the right is one term, so either 
everything converges or everything diverges. 

You can do better than this and use these inequalities to get a quick estimate of the sum of a 
series that would be too tedious to sum by itself. For example 



^ n2 ~ ^ 22 + 32 + ^ n2 

1 4 

This last sum lies between two integrals. 

POO 2 00 ^ /"CXD 2 

73 ^^^>E^> (2.10) 

that is, between 1/3 and 1/4. Now I'll estimate the whole sum by adding the first three terms explicitly 
and taking the arithmetic average of these two bounds. 



00 ^ 

5Z n'- 



1 1 1 /I 1\ , , 

The exact sum is more nearly 1.644934066848226, but if you use brute-force addition of the original 
series to achieve accuracy equivalent to this 1.653 estimation you will need to take about 120 terms. 
This series converges, but not very fast. See also problem 2.24. 

Quicker Comparison Test 

There is another way to handle the comparison test that works very easily and quickly (if it's applicable). 
Look at the terms of the series for large n and see what the approximate behavior of the n^^ term is. 
That provides a comparison series. This is better shown by an example: 

^ — 2n + l/n 
^ 5n^ + sin n 

For large n, the numerator is essentially and the denominator is essentially 5n^, so for large n this 
series is approximately like 



5^2 

More precisely, the ratio of the n^^ term of this approximate series to that of the first series goes to 
one as n — > 00. This comparison series converges, so the first one does too. If one of the two series 
diverges, then the other does too. 

Apply the ratio test to the series for e^. 

e- = yxV^! so ^^ ^'-"'/jk + iy. ^x 
^ ' Uk x^/k\ k + l 
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As A; — > (X) this quotient approaches zero no matter the value of x. This means that the series converges 

for all X. 

Absolute Convergence 

If a series has terms of varying signs, that should help the convergence. A series is absolutely convergent 
if it converges when you replace each term by its absolute value. If it's absolutely convergent then it 
will certainly be convergent when you reinstate the signs. An example of a series that is convergent 
but not absolutely convergent is 

oo ^ 

= ^ " 2 + 3 " • ■ ■ = ^^^^ + = (^-^^^ 

k=l 

Change all the minus signs to plus and the series is divergent. (Use the integral test.) 

Can you rearrange the terms of an infinite series? Sometimes yes and sometimes no. If a series 
is convergent but not absolutely convergent, then each of the two series, the positive terms and the 
negative terms, is separately divergent. In this case you can rearrange the terms of the series to converge 
to anything you want! Take the series above that converges to In 2. I want to rearrange the terms 
so that it converges to \/2. Easy. Just start adding the positive terms until you've passed \/2. Stop 
and now start adding negative ones until you're below that point. Stop and start adding positive terms 
again. Keep going and you can get to any number you want. 

111111 11 
' + 3 + 5-2 + 7 + 9 + n + 13-3"''- 

2.4 Series of Series 

When you have a function whose power series you need, there are sometimes easier ways to the result 
than a straight-forward attack. Not always, but you should look first. If you need the expansion of 
gftx +bx 3bout the origin you can do a lot of derivatives, using the general form of the Taylor expansion. 
Or you can say 

^ax^+bx ^ ^ ^ (q^2 ^ ^ hax^ + hxf + hax'^ + hxf + ■■■ (2.13) 

2 6 

and if you need the individual terms, expand the powers of the binomials and collect like powers of x: 

l + hx + {a + \?ll)x^ + {ah + }?l^)x^ + • • • 

If you're willing to settle for an expansion about another point, complete the square in the exponent 

gOa^+fea; _ ^aix^+bxja) _ ^a{x'^+bx/a+b'^ /4:a?)-b^ /Aa _ ^a{x+b/2a)'^~b^ /Aa _ ^a{x+b/2a)'^ ^-b'^ / ia 

= e-'''/4« [1 + a{x + b/2af + a'^{x + b/2af/2 + • • • ] 

and this is a power series expansion about the point Xq = —b/2a. 

What is the power series expansion of the secant? You can go back to the general formulation 
and differentiate a lot or you can use a combination of two known series, the cosine and the geometric 



series. 



sec a; 



= l+[^x'-^x' + ...] + [^x'-^x^ + ...]'+ ^ ■ ^ 

= l + ^X^+(-iT + (^)V + --- 
1^2 , A 
2!'*' """24 
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This is a geometric series, each of whose terms is itself an infinite series. It still beats plugging into the 
general formula for the Taylor series Eq. (2.5). 
What is 1/sin^x? 



sin^ 



X (x - a;V3! + xV5! ) x^{l - x"^ /V. + /b\ 



1 



3 



= ^ (1 - 3(-xV3! + xV5! -...) + 6(-xV3! + xVs! - . . .f 

_ 1 1 Six 
~ x3 ^ 2a;^ 360 
which is a Frobenius series. 

2.5 Power series, two variables 

The idea of a power series can be extended to more than one variable. One way to develop it is to use 
exactly the same sort of brute-force approach that I used for the one-variable case. Assume that there 
is some sort of infinite series and successively evaluate its terms. 

/(x, y) = A + Bx + Cy + Dx"^ + Exy + Fy"^ + Gx^ + Hx'^y + Ixy"^ + Jy^ ■ ■ ■ 

Include all the possible linear, quadratic, cubic, and higher order combinations. Just as with the single 
variable, evaluate it at the origin, the point (0,0). 

f{0,0) = A + + + --- 

Now differentiate, but this time you have to do it twice, once with respect to x while y is held constant 
and once with respect to y while x is held constant. 

^(x,y) = B + 2Dx + Ey + --- then ^(0,0) = S 
ox ox 

%{x,y) = C + Ex + 2Fy + --- then ^(0,0) = C 
oy oy 

Three more partial derivatives of these two equations gives the next terms. 



dx"^ 
dxdy 

dy"^ 



(x, y) = 2D + QGx + 2Hy ■ ■ 
{x,y) = E + 2Hx + 2ly--- 

(x, y) = 2F + 2lx + 6Jy ■ ■ ■ 



Evaluate these at the origin and you have the values of D, E, and F. Keep going and you have all the 
coefficients. 

This is awfully cumbersome, but mostly because the crude notation that I've used. You can 
make it look less messy simply by choosing a more compact notation. If you do it neatly it's no harder 
to write the series as an expansion about any point, not just the origin. 

oo 

f{x,y)= Aran{x-ar{y-hr (2.15) 

m,n=0 
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Differentiate this m times with respect to x and n times with respect to y, then set x = a and y = b. 
Only one term survives and that is 

Qm+n j 

{a,b) = m\n\Amn 



dx'^dy'' 

I can use subscripts to denote differentiation so that ^ is fx and ^^2^^ is fxxy- Then the 
two-variable Taylor expansion is 

/(x,2/) = /(0)+/,(0)a; + /,(0)y+ 

\ [/xx(0)x2 + 2fxy{^)xy + fyy{^)y^] + 

^ [fxxx{^)X^ + 3fxxy{0)x^y + 3fxyy{0)xy^ + fyyy{0)y^] +■■■ (2.16) 

Again put more order into the notation and rewrite the general form using Amn as 

1 f{m + n)\\ d'^+^'f 
- {m + n)\\ m\n\ J dx^dy^ ^ ^^'^^^ 

That factor in parentheses is variously called the binomial coefficient or a combinatorial factor. Standard 
notations for it are 

m! „ /m\ 

iCn=i (2.18) 



n\{m — n)\ \n ^ 

The binomial series, Eq. (2.4), for the case of a positive integer exponent is 

(1 + x)^ = ^ I I x"^, or more symmetrically 

n=0 

(a + 6)"^ = ^ Pja"6" (2.19) 

n=0 ^'^''^ 

(a + 6)2 = a^ + 2ab + 6^ (a + bf = + 3a% + 3ab'^ + 6^, 
(a + 6)^ = + 4a^6 + Ga^fo^ + 4a53 ^ ^^^^ 

Its relation to combinatorial analysis is that if you ask how many different ways can you choose n 
objects from a collection of m of them, mCn is the answer. 

2.6 Stirling's Approximation 

The Gamma function for positive integers is a factorial. A clever use of infinite series and Gaussian 
integrals provides a useful approximate value for the factorial of large n. 

n\ ~ \/27mn''e"" for large n (2.20) 

Start from the Gamma function of n + 1. 



rco rc 

n\ = r(n + l) = / dtf'e-^ = 

Jo Jo 



dte 



-t+n\nt 



The integrand starts at zero, increases, and drops back down to zero as t — > oo. The graph roughly 
resembles a Gaussian, and I can make this more precise by expanding the exponent around the point 
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where it is a maximum. The largest contribution to the whole integral comes from the region near this 
point. Differentiate the exponent to find the maximum: 



21.06 




d , s n 

( — t + nlntj=— 1 + — = gives t = n 



dt 



Expand about this point 



f(t) = -t + n\nt = f{n) + {t-n)f' {n) + {t - nf f" {n) /2 +••• 

= -n+nlnn + ^ + {t-nf{-n/n^)/2 + ■ ■ ■ 



Keep terms to the second order and the integral is approximately 

^^g-n+nlnn-{t-n)72n ^^ng-n / ^^^-(t-nf/2n ^2.21) 
Jo ^0 

At the lower limit of the integral, at t = 0, this integrand is e^""/^, so if n is even moderately large 
then extending the range of the integral to the whole line — oo to +oo won't change the final answer 
much. 

/oo 
-oo 

where the final integral is just the simplest of the Gaussian integrals in Eq. (1.10). 
To see how good this is, try a few numbers 



n 


n! 


Stirling 


ratio 


difference 


1 


1 


0.922 


0.922 


0.078 


2 


2 


1.919 


0.960 


0.081 


5 


120 


118.019 


0.983 


1.981 


10 


3628800 


3598695.619 


0.992 


30104.381 



You can see that the ratio of the exact to the approximate result is approaching one even though the 
difference is getting very large. This is not a handicap, as there are many circumstances for which this 
is all you need. This derivation assumed that n is large, but notice that the result is not too bad even 
for modest values. The error is less than 2% for n = 5. There are even some applications, especially in 
statistical mechanics, in which you can make a still cruder approximation and drop the factor \j2Tin. 
That is because in that context it is the logarithm of n\ that appears, and the ratio of the logarithms 
of the exact and even this cruder approximate number goes to one for large n. Try it. 

Although I've talked about Stirling's approximation in terms of factorials, it started with the 
Gamma function, so Eq. (2.20) works just as well for Y{n + 1) for any real n: 
r(11.34 = 10.34 + 1) = 8 116 833.918 and Stirling gives 8 051 701. 



2 — Infinite Series 



10 



Asymptotic 

You may have noticed the symbol that I used in Eqs. (2.20) and (2.21). "~" doesn't mean "ap- 
proximately equal to" or "about," because as you see here the difference between nl and the Stirling 
approximation grows with n. That the ratio goes to one is the important point here and it gets this 
special symbol, "asymptotic to." 

Probability Distribution 

In section 1.4 the equation (1.17) describes the distribution of the results when you toss a coin. It's 
straight-forward to derive this from Stirling's formula. In fact it is just as easy to do a version of it for 
which the coin is biased, or more generally, for any case that one of the choices is more likely than the 
other. 

Suppose that the two choices will come up at random with fractions a and b, where a -\-b = 1. 
You can still picture it as a coin toss, but using a very unfair coin. Perhaps a = 1/3 of the time it 
comes up tails and 6 = 2/3 of the time it comes up heads. If you toss two coins, the possibilities are 

TT HT TH HH 

and the fractions of the time that you get each pair are respectively 



a 



2 



ba ab 6^ 



This says that the fraction of the time that you get no heads, one head, or two heads are 

= i/g, 2ab = % 6^ ^ 4/g with total (a + bf = a'^ + 2ab + b"^ = 1 (2.22) 

Generalize this to the case where you throw N coins at a time and determine how often you 
expect to see 0, 1 heads. Equation (2.19) says 



{a + bf = Y,{^^ia%''-^ where 



kj k\{N-k)\ 



When you make a trial in which you toss coins, you expect that the "a" choice will come up N 
times only the fraction of the trials. All tails and no heads. Compare problem 2.27. 

The problem is now to use Stirling's formula to find an approximate result for the terms of this 
series. This is the fraction of the trials in which you turn up k tails and N — k heads. 

k.N-k N\ V2^N^e-^ 



k\{N-ky. fc^e-fc V27r(Ar - k) {N - k)N-ke-(N-k) 

~ V27r^ k{N-k) kk{N-k)N-k 

The complicated parts to manipulate are the factors with all the exponentials of k in them. Pull them 
out from the denominator for separate handling, leaving the square roots behind. 

A;fc(iV-A;)^-W(^-^) 

The next trick is to take a logarithm and to do all the manipulations on it. 



In ^ A; In A; + (A^ - A;) ln(A^ - A;) - A; In a - (A^ - A;) In 6 = f{k) (2.24) 
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The original function is a maximum when this denominator is a minimum. When the numbers N and 
k are big, you can treat k as a continuous variable and differentiate with respect to it. Then set this 
derivative to zero and finally, expand in a power series about that point. 

4rf{k) = lnA; + l-ln(A^-A;)-l-lna + ln6 = 

Ctrh 

In , = In ^ , , = T : k = aN 

N-k 6 N-k b 

This should be no surprise; a is the fraction of the time the first choice occurs, and it says that the 
most likely number of times that it occurs is that fraction times the number of trials. At this point, 
what is the second derivative? 

when fc = aiV, f"(k) = y + ^=^ = -^+ \^ =-^+ ^ ^ 



k N-k aN N-aN aN bN abN 
About this point the power series for f{k) is 

' aN)f'{aN) + ^ 

1 



f{k) = f{aN) + {k- aN)f'{aN) + -{k - aNff"{aN) + 



NlnN+—{k-aN)' + ... (2.25) 

To substitute this back into Eq. (2.23), take its exponential. Then because this will be a fairly sharp 
maximum, only the values of k near to aN will be significant. That allows me to use this central value 
of k in the slowly varying square root coefficient of that equation, and I can also neglect higher order 
terms in the power series expansion there. Let 5 = k — aN. The result is the Gaussian distribution. 



aN{N -aN) N^e^^l^"^^ V2abNTT 

When a = b= 1/2, this reduces to Eq. (1.17). 

When you accumulate N trials at a time (large A^) and then look for the distribution in these 
cumulative results, you will commonly get a Gaussian. This is the central limit theorem, which says 
that whatever set of probabilities that you start with, not just a coin toss, you will get a Gaussian by 
averaging the data. {Not really true. There are some requirements* on the probabilities that aren't 
always met, but if as here the variable has a bounded domain then it's o.k. See problems 17.24 and 
17.25 for a hint of where a naive assumption that all distributions behave the same way that Gaussians 
do can be misleading.) If you listen to the clicks of a counter that records radioactive decays, they sound 
(and are) random, and the time interval between the clicks varies greatly. If you set the electronics to 
click at every tenth count, the result will sound regular, and the time interval between clicks will vary 
only slightly. 



* finite mean and variance 
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2.7 Useful Tricks 

There are a variety of ways to manipulate series, and while some of them are simple they are probably 
not the sort of thing you'd think of until you've seen them once. Example: what is the sum of 

-, 1111 

1 \ 1 ? 

3 5 7 9 

Introduce a parameter that you can manipulate, like the parameter you sometimes introduce to do 
integrals as in Eq. (1.5). Consider the series with the parameter x in it. 

/y*3 

f{x) = x- y + -- y + — - (2.27) 

Differentiate this with respect to x to get 

f(x) = \-x^ +x'^ -x^ ^x^ 



That looks a bit like the geometric series except that it has only even powers and the signs alternate. 
Is that too great an obstacle? As 1/(1 — ,t) has only plus signs, then change x to — ,t, and 1/(1 + x) 
alternates in sign. Instead of a; as a variable, use x^ , then you get exactly what you're looking for. 



f{x) = \-x^ + x^ -x^ + x^ = - 



+ x^ 



Now to get back to the original series, which is /(I) recall, all that I need to do is integrate this 
expression for f'{x). The lower limit is zero, because /(O) = 0. 



/(I) 



dx 



1 + x^ 



tan ^ X 



TT 

4 



This series converges so slowly however that you would never dream of computing tt this way. If you 
take 100 terms, the next term is 1/201 and you can get a better approximation to tt by using 22/7. 
The geometric series is 1 + x + x'^ + x^ + ■ ■ ■, but what if there's an extra factor in front of each 

term? 



fix) 



2 + 3x + 4x'^ + 5x^ + 



Multiply this by x and it is 2x + 3x'^ + ix^ + 5x'^ + • • •, starting to look like a derivative. 

d 



xf{x) = 2x + + 4x'^ + + 



dx 



(x^ + + H ) 



Again, the geometric series pops up, though missing a couple of terms. 



xf{x) 



d^ 
dx 



{l + x + x"^ + x^ -\ 1-x) 



d^ 

dx 



The final result is then 



fix) 



1 

X 



(1-X)2 



1-X 



2-x 

(1^ 



— 1 — X 



ii-xy 



- 1 
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2.8 Diffraction 

When light passes through a very small opening it will be diffracted so that it will spread out in a 
characteristic pattern of higher and lower intensity. The analysis of the result uses many of the tools 
that you've looked at in the first two chapters, so it's worth showing the derivation first. 

The light that is coming from the left side of the figure has a wavelength A and wave number 
k = 27r/A. The light passes through a narrow slit of width = a. The Huygens construction for the 
light that comes through the slit says that you can effectively treat each little part of the slit as if it is 
a source of part of the wave that comes through to the right. (As a historical note, the mathematical 
justification for this procedure didn't come until about 150 years after Huygens proposed it, so if you 
think it isn't obvious why it works, you're right.) 




Call the coordinate along the width of the slit y, where < y < a. I want to find the total light 
wave that passes through the slit and that heads at the angle 6 away from straight ahead. The light 
that passes through between coordinates y and y + dy \s a wave 



Ady cos{kr — cut) 

Its amplitude is proportional to the amplitude of the incoming wave. A, and to the width dy that I am 
considering. The coordinate along the direction of the wave is r. The total wave that will head in this 
direction is the sum (integral) over all these little pieces of the slit. 

Let ro be the distance measured from the bottom of the slit to where the light is received far 
away. Find the value of r by doing a little trigonometry, getting 

r = To — ysinO 

The total wave to be received is now the integral 

sin (A;(ro — y sin 9) — ujt) 



i-a c 
/ Ady cos [k{ro — y sin 9) — out) = A- 
Jq 



-ksm9 



Put in the limits to get 



A 



-ksm9 



[sin [k{ro — asin^) — out) — sin [kro — out)] 



I need a trigonometric identity here, one that you can easily derive with the techniques of complex 
algebra in chapter 3. 



I X — y \ ( X + y 
sin X — sm y = 2 sm ( — - — j cos 



(2.28) 
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Use this and the light amplitude is 

2A . 



; sm 



sm 



ka 



sin j cos [k{rQ — ^ sin 6) — cot 



(2.29) 



The wave is the cosine factor. It is a cosine of {k ■ distance — Lut), and the distance in question 
is the distance to the center of the slit. This is then a wave that appears to be coming from the middle 
of the slit, but with an amplitude that varies strongly with angle. That variation comes from the other 
factors in Eq. (2.29). 

It's the variation with angle that's important. The intensity of the wave, the power per area, 
is proportional to the square of the wave's amplitude. I'm going to ignore all the constant factors, so 
there's no need to worry about the constant of proportionality. The intensity is then (up to a factor) 



sin ((A;a/2)sin6') 



sin^^ 



(2.30) 



For light, the wavelength is about 400 to 700 nm, and the slit may be a millimeter or a tenth of a 
millimeter. The size of ka/2 is then about 

ka/2 = Tca/X ^3-0.1 mm/500 nm ^ 1000 

When you plot this intensity versus angle, the numerator vanishes when the argument of sin^() is mr, 
with n an integer, +, — , or 0. This says that the intensity vanishes in these directions except for = 0. 
In that case the denominator vanishes too, so you have to look closer. For the simpler case that 9 0, 
these angles are 

nil = — smU ^ — U n = ±1, ±2, . . . 
2 2 ' ' 

Because ka is big, you have many values of n before the approximation that sin^^ = 6 becomes invalid. 

You can rewrite this in terms of the wavelength because k = 2n/X. 

27ra^ , 
& = nA/a 



nir 



2X 



or 



What happens at zero? Use power series expansions to evaluate this indeterminate form. The 
first term in the series expansion of the sine is 6 itself, so 

v2 



sin^ ((A;a/2)sin6') 



sm 



20 



{{ka/2)ey 



f ka\' 
\2) 



(2.31) 



What is the behavior of the intensity near 6 = 07 Again, use power series expansions, but keep 
another term 

sine = 6- -6^^ , and {1 + x)" = 1 + ax + ■ ■ ■ 

6 

Remember, ka/2 is big! This means that it makes sense to keep just one term of the sine expansion 
for sinO itself, but you'd better keep an extra term in the expansion of the sm^{ka . . .). 

-\ 2 



sm 



{{ka/2)e) 



1 

1 



ka 



ka , 



ka , 
~2 

1 f ka 



e] + 



ka\- 



6 V 2 



ka . 
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When you use the binomial expansion, put the binomial in the standard form, (1 + x) as in the second 
line of these equations. What is the shape of this function? Forget all the constants, and it looks like 
1 — 9'^. That's a parabola. 



The dots are the points where the intensity goes to zero, nX/a. Between these directions it 
reaches a maximum. How big is it there ? These maxima are about halfway between the points where 
{kasm6)/2 = nvr. This is 



ka 



sm9 = {n+y2)n, n = ±1, ±2, . 



At these angles the value of / is, from Eq. (2.30), 



^ / ka\' 



V 2 y \{2n + l)Ti/2 

The intensity at ^ = is by Eq. (2.31), {ka/2)^, so the maxima off to the side have intensities that 
are smaller than this by factors of 



-^rr- = 0.045, 
97r2/4 



— \— =0.016,... 

257r2/4 



2.9 Checking Results 

When you solve any problem, or at least think that you've solved it, you're not done. You still have to 
check to see whether your result makes any sense. If you are dealing with a problem whose solution is 
in the back of the book then do you think that the author is infallible? If there is no back of the book 
and you're working on something that you would like to publish, do you think that you're infallible? 
Either way you can't simply assume that you've made no mistakes; you have to look at your answer 
skeptically. 

There's a second reason, at least as important, to examine your results: that's where you can 
learn some physics and gain some intuition. Solving a complex problem and getting a complicated 
answer may involve a lot of mathematics but you don't usually gain any physical insight from doing it. 
When you analyze your results you can gain an understanding of how the mathematical symbols are 
related to physical reality. Often an approximate answer to a complicated problem can give you more 
insight than an exact one, especially if the approximate answer is easier to analyze. 

The first tool that you have to use at every opportunity is dimensional analysis. If you are 
computing a length and your result is a velocity then you are wrong. If you have something in your 
result that involves adding a time to an acceleration or an angle to a distance, then you've made a 
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mistake; go back and find it. You can do this sort of analysis everywhere, and it is one technique that 
provides an automatic error finding mechanism. If an equation is dimensionally inconsistent, backtrack 
a few lines and see whether the units are wrong there too. If they are correct then you know that your 
error occurred between those two lines; then further narrow the region where the mistake happened by 
looking for the place at which the dimensions changed from consistent to inconsistent and that's where 
the mistake happened. 

The second tool in your analysis is to examine all the parameters that occur in the result and 
to see what happens when you vary them. Especially see what happens when you push them to an 
extreme value. This is best explained by some examples. Start with some simple mechanics to see the 
procedure. 




Two masses are attached by a string of negligible mass and that is wrapped around a pulley 
of mass M so that it can't slip on the pulley. Analyze them to determine what is wrong with each. 
Assume that there is no friction between mi and the table and that the string does not slip on the 
pulley. 

1112 - mi ,,s m2 , . 1712- M/2 
(a) ax = Q (b) Qx = TTi-Q (c) ax = im-Q 

(a) If mi ^ m2, this is negative, meaning that the motion of mi is being slowed down. But 
there's no friction or other such force to do this. 

OR If mi = m2, this is zero, but there are still unbalanced forces causing these masses to accelerate. 

(b) If the combination of masses is just right, for example mi = 1 kg, m2 = 1 kg, and M = 2 kg, 
the denominator is zero. The expression for ax blows up — a very serious problem. 

OR If M is very large compared to the other masses, the denominator is negative, meaning that is 
negative and the acceleration is a braking. Without friction, this is impossible. 

(c) If M S> mi and m2, the numerator is mostly — Af /2 and the denominator is mostly +M/2. 
This makes the whole expression negative, meaning that mi and m2 are slowing down. There is no 
friction to do this, and all the forces are the direction to cause acceleration toward positive x. 

OR If m2 = M/2, this equals zero, saying that there is no acceleration, but in this system, will 
always be positive. 

The same picture, but with friction //^ between mi and the table. 

, , m2 ... m2 - ^l\,ml , . m2 

a «x = ^ , ..,^ 9 b) ax = ^^^9 (c) dx = ^ r^9 

m2 + Hvmi + M/2 m2 — M/2 m2 + /ikirii — M/2 

(a) If is very large, this approaches zero. Large friction should cause mi to brake to a halt 
quickly with very large negative ax- 

OR If there is no friction, /i.^ = 0, then mi plays no role in this result but if it is big then you know 
that it will decrease the downward acceleration of m2. 

(b) The denominator can vanish. If m2 = M/2 this is nonsense. 

(c) This suffers from both of the difficulties of (a) and (b). 

Trajectory Example 

When you toss an object straight up with an initial speed Vq, you may expect an answer for the motion 
as a function of time to be something like 
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Vy{t) =vo- gt, y{t) = vot - \gt 



(2.32) 



Should you expect this? Not if you remember that there's air resistance. If I claim that the answers are 



Vyit) 



-Vt + [vq + ft e 



-gt/vt] 



(2.33) 



then this claim has to be inspected to see if it makes sense. And I never bothered to tell you what 
the expression "Vt" means anyway. You have to figure that out. Fortunately that's not difficult in this 
case. What happens to these equations for very large time? The exponentials go to zero, so 



-Vt + (f + Wt) ■ 



and 



y 



'Vtt + {Vo + Vt) 



Vt 

9 



Vt is the terminal speed. After a long enough time a falling object will reach a speed for which the force 
by gravity and the force by the air will balance each other and the velocity then remains constant. 
Do they satisfy the initial conditions? Yes: 



Vy{0) = -Vt + {vq +ft)e° 



Vo, y{0) = + {vo + Vt)jil-l) = 



What do these behave like for small time? They ought to reduce to something like the expressions 
in Eq. (2.32), but just as important is to determine what the deviation from that simple form is. Keep 
some extra terms in the series expansion. How many extra terms? If you're not certain, then keep one 
more than you think you will need. After some experience you will usually be able to anticipate what 
to do. Expand the exponential: 



Vy{t) = -Vt + {Vo + Vt) 



1 + 



Z9t 

Vt 



+ 



Vo 



Vo 
Vt 



+ o 1 + 











Vo^ 






Vt 



+ 



+ 



The coefficient of t says that the object is slowing down more rapidly than it would have without air 
resistance. So far, so good. Is the factor right? Not yet clear, so keep going. Did I need to keep terms 
to order t^? Probably not, but there wasn't much algebra involved in doing it, so it was harmless. 
Look at the other equation, for y. 



y{t) = -Vtt + (fo + Vt) 



Vt 

9 



Vot 



qt 1 

— + - 
Vt 2 



1 + 



^^0 



gt 

Vt 

gH^ 



Vt J Vt 



+ 



+ 



Now differentiate this approximate expression for y with respect to time and you get the approximate 
expression for Vy. That means that everything appears internally consistent, and I haven't introduced 
any obvious error in the process of approximation. 

What if the terminal speed is infinite, so there's no air resistance. The work to answer this is 
already done. Expanding e"^*/''* for small time is the same as for large ft, so you need only look back 



2 — Infinite Series 



18 



at the preceding two sets of equations and let — oo. The result is precisely the equations (2.32), 
just as you should expect. 

You can even determine something about the force that I assumed for the air resistance: Fy = 
may = mdvy/dt. Differentiate the approximate expression that you already have for Vy, then at least 



for small t 

d 

^y = '^dt 



1 + — + + — ^ — + • ■ 
/ 2 V vj 

= —m ^1 + — ^ 9 + ■ ■ ■ = —Tng — mgvo/vt + ■ ■ ■ (2.34) 

This says that the force appears to be (1) gravity plus (2) a force proportional to the initial velocity. 
The last fact comes from the factor Vq in the second term of the force equation, and at time zero, that 
is the velocity. Does this imply that I assumed a force acting as Fy = —mg — (a constant times)fy? 
To this approximation that's the best guess. (It happens to be correct.) To verify it though, you would 
have to go back to the original un-approximated equations (2.33) and compute the force from them. 

t>b 



Electrostatics Example 

Still another example, but from electrostatics this time: Two thin circular rings have radii a and b and 
carry charges Qi and Q2 distributed uniformly around them. The rings are positioned in two parallel 
planes a distance c apart and with axes coinciding. The problem is to compute the force of one ring 
on the other, and for the single non-zero component the answer is (perhaps) 



2vr2eo Jo [c^ + {b - af + Aabs\x^ of^' 

Is this plausible? First check the dimensions! The integrand is (dimensionally) l/{c^)^/'^ = l/c^, where 
c is one of the lengths. Combine this with the factors in front of the integral and one of the lengths (c's) 
cancels, leaving Q1Q2/ ^oC^ ■ This is (again dimensionally) the same as Coulomb's law, qiq2/4:Tceor'^ , 
so it passes this test. 

When you've done the dimensional check, start to consider the parameters that control the result. 
The numbers a, b, and c can be anything: small, large, or equal in any combination. For some cases 
you should be able to say what the answer will be, either approximately or exactly, and then check 
whether this complicated expression agrees with your expectation. 

If the rings shrink to zero radius this has a = b = 0, so Fz reduces to 

Q1Q2C r'\Q^ _QiQ2C vr _ Q1Q2 



and this is the correct expression for two point charges a distance c apart. 

If c ^ a and b then this is really not very different from the preceding case, where a and b are 



zero. 



If a = this is 

p Q1Q2C ri^ de _ Q1Q2C Ti/2 _ Q1Q2C 



2vr2eo 7o [c2 + 62]3/2 27r2eo [^2 + 52] 3/2 47reo[c2 + 62] 



3/2 



(2.36) 
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The electric field on the axis of a ring is something that you can compute easily. The only component of 
the electric field at a point on the axis is itself along the axis. You can prove this by assuming that it's 
false. Suppose that there's a lateral component of E and say that it's to the right. Rotate everything 
by 180° about the axis and this component of E will now be pointing in the opposite direction. The 
ring of charge has not changed however, so E must be pointing in the original direction. This supposed 
sideways component is equal to minus itself, and something that's equal to minus itself is zero. 

All the contributions to E except those parallel the axis add to zero. Along the axis each piece 
of charge dq contributes the component 



C 



dq 



The first factor is the magnitude of the field of the point charge at a distance r = Vc^ + and the 
last factor is the cosine of the angle between the axis and r. Add all the dq together and you get Qi. 
Multiply that by Q2 and you have the force on Q2 and it agrees with the expressions Eq. (2.36) 

If c — > then i^^ — > in Eq. (2.35). The rings are concentric and the outer ring doesn't push 
the inner ring either up or down. 

But wait. In this case, where c 0, what if a = 6? Then the force should approach infinity 
instead of zero because the two rings are being pushed into each other. If a = 6 then 

^ Q1Q2C r/' de ^2 

27r2eo Jo [c2 + 4a2 sin^ ^] 
If you simply set c = in this equation you get 

^ Q1Q20 r'^ de 



2vr2eo 7o [4a2 sin^ ^] 



The numerator is zero, but look at the integral. The variable e goes from to 7r/2, and at the end 
near zero the integrand looks like 

1 1 1 



Here I used the first term in the power series expansion of the sine. The integral near the zero end is 
then approximately 

/■■■• de_-i_ 
Jo ~e^~W 

and that's infinite. This way to evaluate Fz is indeterminate: ■ 00 can be anything. It doesn't show 
that this Fz gives the right answer, but it doesn't show that it's wrong either. 

Estimating a tough integral 

Although this is more difficult, even tricky, I'm going to show you how to examine this case for small 
values of c and not for c = 0. The problem is in figuring out how to estimate the integral (2.37) for 
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small c, and the key is to realize that the only place the integrand gets big is in the neighborhood of 
6 = 0. The trick then is to divide the range of integration into two pieces 

n/2 .A .7r/2 



[c2 + 4a2 sin2 6] 



3/2 



For any positive value of A the second piece of the integral will remain finite even as c — > 0. This means 
that in trying to estimate the way that the whole integral approaches infinity I can ignore the second 
part of the integral. Now choose A small enough that for < < A I can use the approximation 
sin 6' = 6, the first term in the series for sine. (Perhaps A = 0.1 or so.) 

de (le , ^ 

for small c, / 1- / — + lower order terms 

■^0 [c2 + 4a2sin2^]^/^ ^ [c^ + 4a2^2]3/2 

This is an elementary integral. Let 6 = (c/2a)tan0. 

de f^' (c/2a)sec2 0ci0 f^' ^ 1 



cos ( 



■ sin A' 



[c2 + 4a2^2]3/2 7o [c2 + c2tan2(/.]3/2 2ac^ Jo ^ 2ac2 

The limit A' comes from A = (c/2a)tanA', so this implies tan A' = 2aA/c. Now given the tangent 
of an angle, I want the sine — that's the first page of chapter one. 

2aA/c 2aA 
sm A - ' - 



7l + (2aA /c)2 Vc2 + 4a2A2 

As c — )■ 0, this approaches one. Put all of this together and you have the behavior of the integral in 
Eq. (2.37) for small c. 

de 1 , , 

1- ^ + lower order terms 

[c2 + 4a2sin2e]'/' 2ac2 

Insert this into Eq. (2.37) to get 

Q1Q2C 1 Q1Q2 



27r2eo 2ac2 An'^eoac 

Now why should I believe this any more than I believed the original integral? When you are very 
close to one of the rings, it will look like a long, straight line charge and the linear charge density on it is 
then A = Qi/lna. What is the electric field of an infinitely long uniform line charge? Er = X/2Tceor. 
So now at the distance c from this line charge you know the i^-field and to get the force on Q2 you 
simply multiply this field by Q2- 

should be -^Q2 = %^Q2 (2.38) 

and that's exactly what I found in the preceding equation. After all these checks I think that I may 
believe the result, and more than that you begin to get an intuitive idea of what the result ought to 
look like. That's at least as valuable. It's what makes the difference between understanding the physics 
underlying a subject and simply learning how to manipulate the mathematics. 
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Exercises 

1 Evaluate by hand cos 0.1 to four places. 

2 In the same way, evaluate tan 0.1 to four places. 

3 Use the first two terms of the binomial expansion to estimate \/2 = y/1 + 1. What is the relative 
error? [(wrong— right)/right] 

4 Same as the preceding exercise, but for Vl-2. 

5 What is the domain of convergence for x — x"^ + — + x^^ — . . . 

oo 

6 Does ^ cos(n) — cos(n + 1) converge? 

n=0 



oo ^ 

7 Does > — = converge? 

n=l ^ 

oo , 

8 Does — ^converge? 



n=l 



2 S 4 

■X> iX-' iXj 



9 What is the domain of convergence for ^ H — - — —r + • • •? 

1-2 2-2^ 3 ■ S'^ 4-4* 

10 From Eq. (2.1), find a series for ttt. 

(1 — xY 

11 If X is positive, sum the series 1 + + e~^^ + e~^^ + • • • 

12 What is the ratio of the exact value of 20! to Stirling's approximation for it? 

13 For the example in Eq. (2.22), what are the approximate values that would be predicted from 
Eq. (2.26)? 

14 Do the algebra to evaluate Eq. (2.25). 

15 Translate this into a question about infinite series and evaluate the two repeating decimal numbers: 
0.444444..., 0.987987987... 

16 What does the integral test tell you about the convergence of the infinite series n^^7 

n What would the power series expansion for the sine look like if you require it to be valid in arbitrary 
units, not just radians? This requires using the constant "C" as in section 1.1. 
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Problems 

2.1 (a) If you borrow $200,000 to buy a house and will pay it back in monthly installments over 30 years 
at an annual interest rate of 6%, what is your monthly payment and what is the total money that you 
have paid (neglecting inflation)? To start, you have payments p with monthly interest i and after 
all payments your unpaid balance must reach zero. The initial loan is L and you pay at the end of 
each month. 

{{L{1 + i) - p){l + i) - p){l + 1) - p ■■■ times =0 

Now carry on and find the general expression for the monthly payment. Also find the total paid. 

(b) Does your general result for arbitrary N reduce to the correct value if you pay everything back at 
the end of one month? + i) = p\ 

(c) For general N, what does your result become if the interest rate is zero? Ans: $1199.10, $431676 

2.2 In the preceding problem, suppose that there is an annual inflation of 2%. Now what is the total 
amount of money that you've paid in constant dollars? That is, one hundred dollars in the year 2010 is 
worth just $100/1. 02^° = $82.03 as expressed in year-2000 dollars. Each payment is paid with dollars 
of gradually decreasing value. Ans: $324211 

2.3 Derive all the power series that you're supposed to memorize, Eq. (2.4). 

2.4 Sketch graphs of the functions 



2.5 The sample series in Eq. (2.7) has a simple graph (x^ between ~L and +L) Sketch graphs of one, 
two, three terms of this series to see if the graph is headed toward the supposed answer. 

2.6 Evaluate this same Fourier series for x^ at x = L; the answer is supposed to be L? . Rearrange 
the result from the series and show that you can use it to evaluate C(2), Eq. (2.6). Ans: 7r^/6 

2.7 Determine the domain of convergence for all the series in Eq. (2.4). 

2.8 Determine the Taylor series for coshx and sinhx. 

2.9 Working strictly by hand, evaluate >70.999. Also \/50- Ans: Here's where a calculator can tell you 
better than I can. 

2.10 Determine the next, x^, term in the series expansion of the secant. Ans: 61x^/720 

2.11 The power series for the tangent is not as neat and simple as for the sine and cosine. You can 
derive it by taking successive derivatives as done in the text or you can use your knowledge of the series 
for the sine and cosine, and the geometric series. 

tanx = = - ^'/^; + • • • = \x- xVS! + • • • Ul + (-xV2! + •••)] 
cosx 1-x2/2!h L ' JL V / n 

Use the expansion for the geometric series to place all the x^, x^, etc. terms into the numerator, 
treating every term after the "1" as a single small thing. Then collect the like powers to obtain the 
series at least through x^. 
Ans: X + + 2x^/15 + 17x7315 + • • • 
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2.12 What is the series expansion for cscx = 1/ sinx? As in the previous problem, use your knowledge 
of the sine series and the geometric series to get this result at least through . Note; the first term 
in this series is l/x. Ans: l/x + x/6 + 7x^/360 + 31xVl5120 H 

2.13 The exact relativistic expression for the kinetic energy of an object with non-zero mass is 

K = mc^['y-l) where 7 = (l - w^/c^)""^^^ 

and c is the speed of light in vacuum. If the speed v is small compared to the speed of light, find 
an approximate expression for K to show that it reduces to the Newtonian expression for the kinetic 
energy, but include the next term in the expansion to determine how large the speed v must be in order 
that this correction term is 10% of the Newtonian expression for the kinetic energy? Ans: v ^ 0.36 c 

2.14 Use series expansions to evaluate 

1 — cosx , ,. s'mkx 
lim — and lim 

x^-ol — cosha; x^o x 

2.15 Evaluate using series; you will need both the sine series and the binomial series. 

lim f — \ \ 

x^o ysm^ X 

Now do it again, setting up the algebra differently and finding an easier (or harder) way. Ans: 1/3 

2.16 For some more practice with series, evaluate 

,• ,2 1 

lim — h 



x^O \X 1 - Vl + X, 

Ans: Check experimentally with a few values of x on a pocket calculator. 

2.17 Expand the integrand to find the power series expansion for 

ln(l +x) = / dt{l + t)-^ 
Jo 

Ans: Eq. (2.4) 

2.18 (a) The error function erf(x) is defined by an integral. Expand the integrand, integrate term 
by term, and develop a power series representation for erf. For what values of x does it converge? 
Evaluate erf(l) from this series and compare it to the result of problem 1.34. (b) Also, as further 
validation of the integral in problem 1.13, do the power series expansion of both sides of the equation 
and verify the expansions of the two sides of the equation agree . 

2.19 Verify that the combinatorial factor mCn is really what results for the coefficients when you 
specialize the binomial series Eq. (2.4) to the case that the exponent is an integer. 

2.20 Determine the double power series representation about (0,0) of l/[(l — x/a)(l ~y/b)~\ 

2.21 Determine the double power series representation about (0,0) of — x/a — y/b) 
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2.22 Use a pocket calculator that can handle 100! and find the ratio of Stirling's approximation to the 
exact value. You may not be able to find the difference of two such large numbers. An improvement 
on the basic Stirling's formula is 

v&n"e-" ( 1 + — 1 
\ 12nJ 

What is the ratio of approximate to exact for n = 1, 2, 10? 
Ans: 0.99898, 0.99948, . . . 

2.23 Evaluate the sum l/n{n + l). To do this, write the single term l/n{n + l) as a combination 
of two fractions with denominator n and {n + 1) respectively, then start to write out the stated infinite 
series to a few terms to see the pattern. When you do this you may be tempted to separate it into two 
series, of positive and of negative terms. Examine the problem of convergence and explain why this is 
wrong. Ans: 1 

2.24 (a) You can sometimes use the result of the previous problem to improve the convergence of a 
slow-converging series. The sum l/'^^ converges, but not very fast. If you add zero to it you don't 
change the answer, but if you're clever about how you add it you can change this into a much faster 
converging series. Add 1 — + 1) to this series and combine the sums, (b) After Eq. (2.11) it 
says that it takes 120 terms to get the stated accuracy. Verify this. For the same accuracy, how many 
terms does this improved sum take? Ans: about 8 terms 

2.25 The electric potential from one point charge is kq/r. For two point charges, you add the 
potentials of each: kqi/ri + kq2/r2- Place a charge —q at the origin; place a charge +q at position 
{x,y,z) = (0,0, a). Write the total potential from these at an arbitrary position P with coordinates 
(x, y, z). Now suppose that a is small compared to the distance of P to the origin (r = \/ x"^ + y"^ + z'^^ 
and expand your result to the first non-vanishing power of a, or really of a/r. This is the potential 
of an electric dipole. Also express your answer in spherical coordinates. See section 8.8 if you need. 
Ans: kqa cosO/r"^ 

2.26 Do the previous problem, but with charge —2q at the origin and charges +q at each of the two 
points (0,0, a) and (0,0, —a). Again, you are looking for the potential at a point far away from the 
charges, and up to the lowest non-vanishing power of a. In effect you're doing a series expansion in 
a/r and keeping the first surviving term. Also express the result in spherical coordinates. The angular 
dependence should be proportional to P2{cos6) = — ^, a "Legendre polynomial." The r 
dependence will have a l/r^ in it. This potential is that of a linear quadrupole. 

2.27 The combinatorial factor Eq. (2.18) is supposed to be the number of different ways of choosing 
n objects out of a set of m objects. Explicitly verify that this gives the correct number of ways for 
m = 1, 2, 3, 4. and all n from zero to m. 

2.28 Pascal's triangle is a visual way to compute the values of mCn- Start with the single digit 1 on 
the top line. Every new line is computed by adding the two neighboring digits on the line above. (At 
the end of the line, treat the empty space as a zero.) 

1 

1 1 
1 2 1 



13 3 1 
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Write the next couple of lines of the triangle and then prove that this algorithm works, that is that the 
m*^ row is the rnCn, where the top row has m = 0. Mathematical induction is the technique that I 
recommend. 

2.29 Sum the series and show 

12 3 

2! ^ 3! ^ 4! ^ ■ ■ ■ ~ ^ 

2.30 You know the power series representation for the exponential function, but now apply it in a 
slightly different context. Write out the power series for the exponential, but with an argument that 
is a differential operator. The letter h represents some fixed number; interpret the square of d/dx as 
d^/dx"^ and find 

e'^^.fix) 

Interpret the terms of the series and show that the value of this is f{x + h). 

2.31 The Doppler effect for sound with a moving source and for a moving observer have different 
formulas. The Doppler effect for light, including relativistic effects is different still. Show that for low 
speeds they are all about the same. 

The symbols have various meanings: v is the speed of sound in the first two, with the other terms 
being the velocity of the observer and the velocity of the source. In the third equation c is the speed 
of light and v is the velocity of the observer. And no, 1 = 1 isn't good enough; you should get these 
at least to first order in the speed. 

2.32 In the equation (2.30) for the light diffracted through a narrow slit, the width of the central 
maximum is dictated by the angle at the first dark region. How does this angle vary as you vary the 
width of the slit, a? What is this angle if a = 0.1 mm and A = 700 nm? And how wide will the central 
peak be on a wall 5 meters from the slit? Take this width to be the distance between the first dark 
regions on either side of the center. 

2.33 An object is a distance d below the surface of a medium with index of refraction n. (For example, 
water.) When viewed from directly above the object in air (i.e. use small angle approximation), the 
object appears to be a distance below the surface given by (maybe) one of the following expressions. 
Show why most of these expressions are implausible; that is, give reasons for eliminating the wrong 
ones without solving the problem explicitly. 

{I) d^/l + ri^/n (2) dn/\/l + n^ (3) nd (4) d/n (5) dn^ / \/ 1 + 




a. 



y 



□ 



2.34 A mass mi hangs from a string that is wrapped around a pulley of mass M. As the mass mi 
falls with acceleration ay, the pulley rotates. An anonymous source claims that the acceleration of mi 
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is one of the following answers. Examine them to determine if any is plausible. That is, examine each 
and show why it could not be correct. NOTE: solving the problem and then seeing if any of these agree 
is not what this is about. 

(1) Qy = Mg/{mi - M) (2) ay = Mg/{mi + M) (3) ay = mig/M 

2.35 Light travels from a point on the left (p) to a point on the 
right {q), and on the left it is in vacuum while on the right of 
the spherical surface it is in glass with an index of refraction n. 
The radius of the spherical surface is R and you can parametrize 
the point on the surface by the angle 9 from the center of the 
sphere. Compute the time it takes light to travel on the indicated 
path (two straight line segments) as a function of the angle 9. 
Expand the time through second order in a power series in 9 and 
show that the function T{9) has a minimum if the distance q is 

small enough, but that it switches to a maximum when q exceeds a particular value. This position is 
the focus. 

2.36 Combine two other series to get the power series in 9 for In(cos^). 
Ans:-i^2-^^4-^^6 + ... 

2.37 Subtract the series for ln(l — x) and ln(l + x). For what range of x does this series converge? 
For what range of arguments of the logarithm does it converge? 

Ans: — 1 < a; < 1, < arg < oo 




2.38 A function is defined by the integral 



dt 



1-^2 



Expand the integrand with the binomial expansion and derive the power (Taylor) series representation 
for / about x = 0. Also make a hyperbolic substitution to evaluate it in closed form. 

2.39 Light travels from a point on the right {p), hits a spherically 
shaped mirror and goes to a point {q). The radius of the spherical 
surface is R and you can parametrize the point on the surface by 
the angle 9 from the center of the sphere. Compute the time 
it takes light to travel on the indicated path (two straight line 
segments) as a function of the angle 9. 

Expand the time through second order in a power series in 9 and 
show that the function T{9) has a minimum if the distance q is 
small enough, but that it switches to a maximum when q exceeds 
a particular value. This is the focus. 

2.40 (a) The quadratic equation ax'^ + 6a; + c = is almost a linear equation if a is small enough: 
hx + c = {)^x = —cjh. You can get a more accurate solution iteratively by rewriting the equation as 




^ b 

Solve this by neglecting the second term, then with this approximate Xi get an improved value of the 
root by 

c an 

^'=-b-r' 
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and you can repeat the process. For comparison take the exact solution and do a power series expansion 
on it for small a. See if the results agree. 

(b) Where does the other root come from? That value of x is very large, so the first two terms in 
the quadratic are the big ones and must nearly cancel, ax^ + bx = so x = —b/a. Rearrange the 
equation so that you can iterate it, and compare the iterated solution to the series expansion of the 
exact solution. 

b c 
a ax 

Solve O.OOlx^ + x + 1 = 0. Ans: Solve it exactly and compare. 

2.41 Evaluate the limits 

, , sinx — tanx . ,. sinx — tanx , , ,. sinx — tanx 
(a) lim , (b) lim k , (c) lim 5 

a:-s>0 X x^O X^ x^O X'^ 

Ans: Check with a pocket calculator for x = 1.0, 0.1, 0.01 

2.42 Fill in the missing steps in the derivation of Eq. (2.26). 

2.43 Is the result in Eq. (2.26) normalized properly? What is its integral d6 over all 57 Ans: 1 

2.44 A political survey asks 1500 people randomly selected from the entire country whom they will 
vote for as dog-catcher-in-chief. The results are 49.0% for T.I. Hulk and 51.0% for T.A. Spiderman. 
Assume that these numbers are representative, an unbiased sample of the electorate. The number 
0.49 X 1500 = aN is now your best estimate for the number of votes Mr. Hulk will get in a sample 
of 1500. Given this estimate, what is the probability that Mr. Hulk will win the final vote anyway? 
(a) Use Eq. (2.26) to represent this estimate of the probability of his getting various possible outcomes, 
where the center of the distribution is at A; = aN. Using 6 = k — aN, this probability function is 
proportional to exp ( — 6^/2abN^ , and the probability of winning is the sum of all the probabilities of 
having k > N/2, that is, C') What would the answer be if the survey had asked 150 or 15000 

people with the same 49-51 results? Ans: (a) ^ [l - erf {^N/2ab{l - a))] . 22%, (b) 40%, 0.7% 

2.45 For the function defined in problem 2.38, what is its behavior near x = 1? Compare this result 
to equation (1.4). Note: the integral is + J^. Also, I - t"^ = {I + t){l - t), and this « 2(1 - t) 
near 1. 

2.46 (a) What is the expansion of 1/(1 -|- t^) in powers of t for small t. (b) That was easy, now what 
is it for large t? In each case, what is the domain of convergence? 

2.47 The "average" of two numbers a and b commonly means {a + b)/2, the arithmetic mean. There 
are many other averages however, (a, 6 > 0) 

M,(a,6)= [(a'^ + 6")/2]'/" 

is the n^^ mean, also called the power mean, and it includes many others as special cases, n = 2: 
root-mean-square, n = —1: harmonic mean. Show that this includes the geometric mean too: Vab = 
liuin-^Q Mn{a,b). It can be shown that dMn/dn > 0; what inequalities does this imply for various 
means? Ans: harmonic < geometric < arithmetic < rms 

2.48 Using the definition in the preceding problem, show that dMn/dn > 0. [Tough!] 
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2.49 In problem 2.18 you found the power series expansion for the error function — good for smal 
arguments. Now what about large arguments? 



erf(a;) = ^ / dte"*' = 1 - ^ / dte"*' = 1 - ^ / ^ -te 



X 



Notice that you can integrate the te * factor explicitly, so integrate by parts. Then do it again and 
again. This provides a series in inverse powers that allows you evaluate the error function for large 
arguments. What is erf(3)? Ans: 0.9999779095 See Abramowitz and Stegun: 7.1.23. 

2.50 A friend of mine got a different result for Eq. (2.35). Instead of s\v? 6 in the denominator, he 
got a sin^. Analyze his answer for plausibility. 

2.51 Find the minimum of the function /(r) = ar+&/r for a, 6, r > 0. Then find the series expansion 
of / about that point, at least as far as the first non-constant term. 

2.52 In problem 2.15 you found the limit of a function as x — )• 0. Now find the behavior of the same 
function as a series expansion for small x, through terms in . Ans: | + j^^x'^. To test whether this 
answer or yours or neither is likely to be correct, evaluate the exact and approximate values of this for 
moderately small x on a pocket calculator. 

2.53 Following Eq. (2.34) the tentative conclusion was that the force assumed for the air resistance was 
a constant times the velocity. Go back to the exact equations (2.33) and compute this force without 
approximation, showing that it is in fact a constant times the velocity. And of course find the constant. 

2.54 An object is thrown straight up with speed Vq. There is air resistance and the resulting equation 
for the velocity is claimed to be (only while it's going up) 

„ m _„ Vo-Vtian{gt/v^) 

Vy[T) — Vt- 



Vt + Vq tan{gt/vt) 

where Vt is the terminal speed of the object after it turns around and has then been falling long enough, 
(a) Check whether this equation is plausible by determining if it reduces to the correct result if there 
is no air resistance and the terminal speed goes to infinity, (b) Now, what is the velocity for small 
time and then use Fy = may to infer the probable speed dependence of what I assumed for the air 
resistance in deriving this expression. See problem 2.11 for the tangent series, (c) Use the exact Vy{t) 
to show that no matter how large the initial speed is, it will stop in no more than some maximum time. 
For a bullet that has a terminal speed of 100 m/s, this is about 16s. 

2.55 Under the same circumstances as problem 2.54, the equation for position versus time is 

-^i^f VtCos{gt/vt) +Vosm{gt/vt] 



9 \ 

(a) What is the behavior of this for small time? Analyze and interpret what it says and whether it 
behaves as it should, (b) At the time that it reaches its maximum height {Vy = 0), what is its position? 
Note that you don't need to have an explicit value of t for which this happens; you use the equation 
that t satisfies. 

2.56 You can get the individual terms in the series Eq. (2.13) another way: multiply the two series: 

^ax'^+bx _ ^ax^ ^bx 



Do so and compare it to the few terms found after (2.13). 



